Low temperature relaxational dynamics of the Ising chain in a transverse field 
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We present asymptotically exact results for the real time order parameter correlations of a class 
of d = 1 Ising models in a transverse field at low temperatures (T) on both sides of the quantum 
critical point. The correlations are a product of a T-independent factor determined by quantum 
effects, and a T-dependent relaxation function which comes from a classical theory. We confirm 
our predictions by a no-free-parameter comparison with numerical studies on the nearest neighbor 
spin-1/2 model. 



Real time, non-zero temperature (T), correlation func- 
tions of quantum many body systems are directly related 
to the observables of many experiments. For strongly in- 
teracting systems, there are few quantitative results on 
the relaxation and transport processes that are believed 
to occur at long times at any T ^ 0. Monte Carlo and 
perturbative methods work best in imaginary time, but 
the analytic continuation to real time is most dangerous, 
and often fails, in the low frequency limit. At special con- 
formally invariant points in dimension d = 1, real time 
correlations describing relaxation of an order parameter 
can be computed exactly. Among systems in arbitrary 
d, which are tuned across a quantum critical point by 
a variable coupling, there is only one for which reliable 
results are available: the d = 1 impenetrable Bose gas, 
whose correlators were determined by ^ profound and 
sophisticated inverse scattering analysis^. 

In this paper, we study real time, T / correlations 
of the d = 1 Ising model in a transverse field. We use 
a novel semiclassical method to obtain the exact asymp- 
totics of order parameter correlations in the two low T 
regions on either side of its quantum critical point. Our 
main new result, given in Eqs. (^), (||) below, is that in 
these low T regimes, the spin correlation function can be 
expressed as a product of two factors. One, arising from 
quantum effects, gives the T = value of the correlation 
function, and the other, which comes from a classical 
theory, describes the effects of temperature. Combined 
with earlier resultsotj, our new results give a complete de- 
scription of time-dependent correlations in all the distinct 
limiting regions of the Ising model, and exhibit, simply 
and clearly, the crossovers in the roles of quantum and 
thermal fluctuations in the relaxational dynamics. 

For the Ising chain with only nearest neighbojjj-jex- 
change, we also use the free fernaipn representationau to 
obtain accurate, numerical dataaQ for real time spin cor- 
relations for systems with up to 512 spins: these results 
are in excellent agreement with the asymptotic theoreti- 
cal results even at relatively short times and distances. 



We also use our semiclassical method to obtain corre- 
lators of the d = 1 impenetrable Bose gas in a certain low 
T regime; our results here are in agreement with earlier 
workLH, although our approach is much simpler and more 
physically transparent. Our method should also apply to 
other d = 1 quantum models (like the non-linear sigma 
or sine-Gordon) with an excitation gap. 

We consider the Hamiltonian 



(1) 



where of, af are Pauli matrices on a chain of sites i, Ji 
(> 0) are short-ranged exchange constants, and g (> 0) is 
the transverse field. The ground state of Hi is expected 
to have long range order with Nq = (a z ) ^ for all 
g <Qc, and a gap to all excitations for g ^ gja, see 
Fig [n. Away from g = g c , the low lying states consist 
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FIG. 1. Finite T phase diagram of Hi. The dashed lines 
are crossover boundaries at T ~ |mc 2 |. 

of a stable particle with energy e p at momentum p, and 
n > 1 particle continua at higher energy. We also have 
E-p = e p , and expect that e p has a minimum at p = 0. 
Consequently, for smallp we parameterize e p — ((mc 2 ) 2 + 
c 2 p 2 + 0(p 4 )) 1 ^ 2 , which defines the 'velocity' c > and 
the 'mass' m; we choose m > (m < 0) for g < g c 
(g > g c ). The critical point at g = g c is described by 
a continuum theory with dynamic exponent z = 1, and 
correlation length exponent v = 1; hence c ~ constant 
for g near g c , while m vanishes linearly, |m| ~ \g — g c \, 
with possibly a different slope on the two sides. 
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The model with Je>i = is integrableH, and all pa- 
rameters are known exactly: in units where Ti = ks = 1 
we have g c = J±, No = (1 — (g/Ji) 2 ) 1 ^ and e p = 
2(Jf + g 2 — 2Jigcos{pa)) 1 / 2 , where a is the lattice spac- 
ing, so mc 2 = 2(Ji — g). However, the existence of these 
stable particles is not a special feature of this integrable 
point. Indeed, for g <C g c , a particle has the simple 
weak-coupling interpretation as the boundary between 
domains with opposite orientations in (cr z ). Conversely, 
for g 3> <?e, the ground states has all the spins oriented in 
the +x direction, and the particle is a — x spin hopping 
from site to site. We expect that these interpretations 
remain correct at the lowest energies, as we approach g c 
from either side, but not at g = g c . 

In the above interpretations, the particles are evidently 
bosons. They have short-range repulsive interactions, 
which, by an elementary calculation! implies that the 
two-particle S'-matrix S pp ' approaches —1 for pa, p 1 a <C 1. 
The integrable model has S pp i = — 1 for all p, p' ', and this 
is often used to obtain a free-fermion description after a 
non-local gauge transformation: we shall not use this 
transformation in our analytical computations. 

Our new analytical results are in the two low T regimes 
on either side of g c , with T <C |mc 2 |. In these regimesta, 
the density, p, of thermally excited particles ~ e _ ' mc '/ T , 
is exponentially small, their mean spacing is much larger 
than their thermal de Broglie wavelength, ~ (2mT) -1 / 2 . 
These are the standard conditions for classical behavior 
in which p — J dp/ {2it)e~^t T . 

We now derive results in the "renormalized classical" 
(RC) region g < g c , T <C mc 2 . Consider the correlation 
function 

C( Xi ,t) = Tr (e- Hl l T e iHlt ale- iHlt a{\ /Z, (2) 

where Xi — ia, Z = Tre~ Hl / T , This can be evaluated 
in terms of the trajectories of the dilute gas of classical 
particles noting that, since the particles physically repre- 
sent domain walls, a z changes sign every time a particle 
goes by. Observe that the classical trajectories remain 
straight lines across collisions because the momenta be- 
fore and after the collision arc the same in d — 1. This im- 
plies that the trajectories are simply independently dis- 
tributed straight lines, placed with a uniform density p 
along the x axis, with an inverse slope v p = de p /dp, and 
with their momenta chosen with the Boltzmann proba- 
bility density e~ e "/ T / p (Fig 0). This heuristic argument, 
based on the classical picture, can be justified by tak- 
ing the semi-classical (stationary phase) limit of the the 
double-time ('Keldysh') path integral], in which each col- 
lision appears in both the forward and backward paths 
(generated by the e~ lHlt and the e lHlt in (0J respec- 
tively) and therefore contributes the factor | S'pp' | 2 = 1. 

Computing C(x, t) is now an exercise in classical prob- 
abilities. The value of a z (0, 0)a z (x, t) is the square of 
the magnetization renormalized by quantum fluctuations 
(Nq), times (—1) if the number of trajectories intersect- 
ing the dashed line in Fig. is odd. Consider an Ising 
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FIG. 2. A typical semiclassical contribution to the double 
time path integral for C'(x,t). Full lines are thermally excited 
particles which propagate forward and backward in time. The 
± signs are significant only for g < g c and denote the orienta- 
tion of the order parameter. For g > g c , the dashed line is a 
particle propagating only forward in time from (0, 0) to (x,t). 

system of size L 3> and let it contain N thermally 
excited particles; then p = N/L. Let the probability that 
any given trajectory intersect the dashed line = q; then 
the probability only a given set of k lines will intersect is 
= q k {l — q) N ~ k . Summing over all possibilities, we have 

N 

C(x,t) = iVo ^2(-l) k q k 0- - q) N ~ k N\/(k\(N - k)\) 

k=0 

= N 2 (l-2q) N ^N 2 e- 2 ^ N . (3) 

The last step holds because, as we shall now compute, 
q <C 1. First, for t — 0, as the density of trajectories 
along the t = axis is uniform, q = \x\/L. For t ^ 0, 
consider first all trajectories with a fixed momentum p: 
they will intersect the dashed line if their intersection 
with t = axis is between ordinates and x — v p t, and 
so q = \x — v p t\/L. Averaging over ^11 p, and inserting in 
(||) , we get one of our main resultstl: 

C(x, t) = N 2 R(x, t) (RC region), (4) 

where R(x, t) = exp (- J — e" £p/T \x - v p t\) . (5) 

Notice that R(x,0) = e — I ^ I an d i?(0,<) = e -^/ T , but 
the general behavior is more complicated. The corre- 
lation length £ = X/(2p). Remarkably, we find from 
(H) that the correlation time, r, is independent of the 
functional form of e p and depends only on the gap: 
t = {n/2T)e mc2 / J '. In the T -> scaling limit, in which 
£ = (7r/2mT) 1 / 2 e mc2 / T , R(x,t) obeys the scaling form 
R(x,t) = <f)(\x\/£, \t\/r), where the scaling function <f> is 
given by 

In (j){x, t) = -x erf (j^J ~le-^'^\ (6) 

Similar classical scaling forms have been discussed ear- 
lierO, but it was incorrectly conjectured that the SGai- 
ing functions would be those of the Glauber modcM; 
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FIG. 3. The points show numerical data for the nearest 
neighbor model, in the RC region, obtained for a lattice size 
L = 256 with free boundary conditions. This is compared 
with the theoretical prediction in Eqs. (Q) and (|^), in which 
the full lattice dispersion, e p was used. In the numerics, the 
two sites were chosen to be as close as possible to the center 
of the lattice. 



Glauber dynamics does not conserve total energy and 
momentum, and these conservation laws have played a 
crucial role in the kinematic constraints on particle col- 
lisions. 

We have also investigated the nearest neighbor model 
numerically by the mapping to free fermionsutj. We gen- 
eralized an earlier .studytj of equal time properties to dy- 
namical quantities^ for the case of free boundary condi- 
tions; details will be presented elsewhere. We find that 
in the RC regime the imaginary part of C(x,t) is much 
smaller than its real part, which suggests that the dy- 
namics is indeed classical as argued above. In Fig. || 
we show data for g = 0.6, T = 0.3 and x — 20 for a 
lattice size L = 256. We took J\ = 1 so g c = 1 and 
mc 2 — 2(1 — g) = 0.8. For comparison we show the the- 
oretical prediction from Eqs. (|j) and (||), in which we 
used the full lattice dispersion relation. The agreement 
is remarkably good. It is interesting to note that because 
there is a maximum velocity, v max , (on the lattice this 
is given by v max = 2J x g for g < g c and v max = 2J X for 
g > g c , whereas in the continuum v max = c), Eq. (j^) 
predicts that C(x,t") should be independent of t for 
\t\ < \x\jvmax, and the numerical results show this very 
clearly. The inset, with a much increased vertical scale, 
gives an idea of how small are the deviations between the 
theory and numerics. We find that the agreement is also 
excellent even at r = 0, for t > 1. 

Now we turn to the "quantum disordered" (QD) region 
g > g c , m < 0, T <C |toc 2 |. The operator a z flips spins 
between the ±x directions, and the large g picture, noted 
earlier, then suggests that <j z is the sum of a creation and 
annihilation operator for the particles. As a result, the 



T = the spectral density, obtained from the Fourier 
transform of C(x,t), has a contribution ~ S(u — e p ) asso- 
ciated with the stable particle. Higher order corrections 
in 1/g, or the form factor expansiontll on the continuum 
theory valid close to <? c , show that the next contribution 
to the spectral density is a continuum above the 3 par- 
ticle threshold. Here we will focus exclusively on how 
the one particle pole broadens as T becomes non-zero. 
We define K(x,t) = C(x,t) at T = 0, and dropping the 
multi-particle terms, we have 



K(x,t) 



1 1 D ^ 



(7) 



where D(p) is a form factor. For the general lattice 
model, D(jp) is not known; neglecting the multi-particle 
terms, D(p) is seen from (0) to be the spatial Fourier 
transform of C(x, 0) i.e. the structure factor. For the 
continuum theory, we have D(p) — Ac/ (2^ where A is 
the dimcnsionlcss quasiparticle amplitudetil (for an un- 
derlying integrable lattice model A = 2(\mc 2 \/ J1) 1 / 4 ), 
whence K(x,t) = AK (mc(x 2 - c 2 i 2 ) 1 / 2 )/(27r), with K 
the modified Bessel function. 

Now we consider T ^ in the semiclassical approxima- 
tion. A typical set of paths contributing to the Keldysh 
path integral is still given by Fig ||, but its physical in- 
terpretation is now very different. The dashed line now 
represents the trajectory of a particle created at (0, 0) 
and annihilated at (x,t), and ± signs in the domains 
should be ignored. In the absence of any other particles 
this dashed line would contribute K(x, t) to C(x, t). The 
scattering off the background particles (the full lines in 
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introduces factors of the S'-matrix element S v 



le dashed line only propagates forward in time, the 
S'-matrix elements for collisions between the dashed and 
full lines (and only these) are not neutralized by a com- 
plex conjugate partner. Using the low momentum value 
Sp P > = — 1, we see that the contribution to C(x, t) equals 
(— l) ne K(x, t) where 111 is the number of full lines inter- 
secting the dashed line. The (— 1)™ £ is precisely the term 
that appeared in the RC region, although for very dif- 
ferent reasons. We can carry out the averaging over all 
trajectories as before, and obtain our main result 



C{x, t) = K(x, t)R(x, t) (QD region) 



(8) 



where K(x, t) is given by Eq. (^) and R{x, t) by Eq. (||). 
While (0) is valid for all x,t, relaxation due to classical 
particles makes sense only within the 'light cone', and so 
strictly speaking, (@) requires ct 3> x. Outside the light 
cone, K(x,t) decays exponentially to zero on the short 
length scale ~ 1/mc at which R ss 1 and T-dependent 
effects are not expected to be large: so it is reasonable 
to use (||) except, perhaps, for x,t extremely small. 

The result (g) clearly displays the separation in scales 
at which quantum and thermal effects act. Quantum 
fluctuations determine the oscillatory, complex function 
K(x,t), which gives the T — value of C(x,t). Expo- 
nential relaxation of spin correlations occurs at longer 
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FIG. 4. The numerical data for the nearest neighbor 
model with Ji = 1, in the QD region, obtained for a lattice 
size L — 512 with free boundary conditions. Also shown is the 
theoretical prediction from Eqs. (^), (JE]) and (Q), with D(p) 
determined in the continuum theory. The numerical data has 
a "ringing" at high frequency, which we believe is due to the 
upper energy cut-off in the dispersion relation, that is present 
in the lattice model but not in the continuum model used 
for the theory. The envelope of the numerical curve fits the 
theoretical prediction well (but not perfectly on the scale of 
resolution of the figure). The insert shows a part of the same 
data with an enlarged scale for clarity. 



of the g = g c critical FWHpt. An explicit expression for 
C(x, t) was given earlierEm, and has exponential decay at 
a single time scale ~ 1/T in both its real and imaginary 
parts. There is no clear separation between the contri- 
butions of thermal and quantum fluctuations since both 
are effective at this scale. Hence there-can be no effective 
classical model describing relaxationlia. For complete- 
ness, we also note the non-universal lattice high T region 
T J\,g\ time-dependent correlators at T = oo were 
obtained earlieia and show a Gaussian decay in time. 

Before concluding, we emphasize that although we 
compared the analytical results with numerical data on 
an integrable model, the main result does not depend on 
integrability. The simplification of the integrable model 
is that the parameters in the theory, m, c, A and Nq, are 
known exactly, whereas in general, they would be phe- 
nomenological parameters of a low energy theory. The 
equality of the classical relaxation functions in the RC 
and QD regions is surely related to the self-duality of the 
Ising critical theory, and is a special feature of this model. 
Finally, it is interesting to speculate that Eq. (||) may also 
be true for d > 1 because of the separation of time scales 
between the quantum and thermal fluctuations. 

We thank S. Majumdar, K. Damle, and T. Senthil for 
helpful discussions. This research was supported by NSF 
Grants No DMR 96-23181 and DMR 94-11964. 



scales ~ £, r, and is controlled by the classical motion of 
particles and a purely real relaxation function R(x,t). 

In Fig. ^ we compare the predictions of Eq. (||) with 
numerical results on a lattice of size L = 512. The the- 
oretical curve was determined from the continuum ex- 
pression for K(x,t), but the full lattice form for e p was 
used in Eq. (||) to determine R(x,t). With the parame- 
ters used, \mc 2 \ = 0.2, which is greater than T (= 0.1) 
as required to be in the QD region. The theory predicts 
that lmC(x 1 t) = for \t\ < \x\/v max , with a singularity 
in both real and imaginary parts at \t\ — \x\/v max . At 
longer time both parts oscillate and decay. Part of this 
decay comes from the \j \ft decay of the Bessel functions, 
but for \t\ > t this is dominated by the exponential decay 
from R(x,t). The theory agrees well with the numerics; 
some differences are visible for small x, outside the light 
cone, but this is outside the domain of validity of (@). 

Precisely the same semiclassical arguments can also 
be applied to the d — 1 impenetrable gas of bosons of 
mass tub, in a chemical potential /i < 0, with T <C 
This system has S- matrix elements S pp ' = —1, and its 
single particle propagator will be given by (||J|,@), with 
e p = —fi + p 2 /2mB and D(p) = 1. This propagator has 
also been computed by the inverse scattering methodEI, 
and the T <C limit of their result agrees with (||). 

Turning to the "quantum critical" region \mc 2 \ «T« 
Ji,<?, the dynamics is now given by finite T correlators 
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